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ABSTRACT 

In this paper we develop the jet model of Potter & Cotter (2012) to include a magnetically 
dominated accelerating parabolic base transitioning to a slowly decelerating conical jet with a 
geometry set by recent radio observations of M87. We conserve relativistic energy-momentum 
and particle number along the jet and calculate the observed synchrotron emission from the jet 
by calculating the integrated line of sight synchrotron opacity through the jet in the rest frame 
of each section of plasma. We calculate the inverse-Compton emission from synchrotron, 
CMB, accretion disc, starlight, broad line region, dusty torus and narrow line region photons 
by transforming into the rest frame of the plasma along the jet. 

We fit our model to simultaneous multi-wavelength observations of the Compton- 
dominant FSRQ type blazar PKS0227-369, with a jet geometry set by M87 and an accelerating 
bulk Lorentz factor consistent with simulations and theory. We investigate models in which the 
jet comes into equipartition at different distances along the jet and equipartition is maintained 
via the conversion of jet bulk kinetic energy into particle acceleration. We find that the jet 
must still be magnetically dominated within the BLR and cannot be in equipartition due to the 
severe radiative energy losses. The model fits the observations, including radio data, very well 
if the jet comes into equipartition outside the BLR within the dusty torus (1.5pc) or at further 
distances (34pc). The fits require a high power jet with a large bulk Lorentz factor observed 
close to the line of sight, consistent with our expectations for a Compton-dominant blazar. We 
find that our fit in which the jet comes into equipartition furthest along the jet, which has a 
jet with the geometry of M87 scaled linearly with black hole mass, has an inferred black hole 
mass close to previous estimates. This implies that the jet of PKS0227 might be well described 
by the same jet geometry as M87. 

Key words: Galaxies: jets, galaxies: active, radiation mechanisms: non-thermal, radio contin- 
uum: galaxies, gamma-rays: galaxies. 



1 INTRODUCTION 

Blazars are thought to be active galactic nuclei (AGN) with rela- 
tivistic jets close to our line of sight. The non-thermal synchrotron 
and inverse-Compton emission from the jet are Doppler-boosted so 
that blazars are the most luminous and variable of AGN, with appar- 
ent superluminal motion visible in their jets. The boosted jet emis- 
sion dominates the spectrum of blazars with synchrotron emission 
extending from radio to UV/x-rays and inverse-Compton emission 
from x-rays to very high energy 7-rays. Blazars are observed to 
produce radio emission which is flat or nearly flat in flux. 

Radio observations show that AGN jets are approximately 
continuous axisymmetric plasma jets. Recent observations show 
that on large scales jets are approximately conical in geometry and 
measurements of the frequency dependent core-shift of M87 and 
other AGN favour a conical jet with conserved magnetic energy. 
VLBI observations of M87 by |Hada et al] ( |2011| > show the blunt 
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base of the jet is within 14 — 23Rs of the central black hole. Ob- 
servations of the geometry of the jet in M87 by Asada & Na kamura] 
( [201 1| | over a large range in length scales show that the jet is in fact 
parabolic close to the base and becomes conical at 10 5 i?s. These 
observations are supported qualitatively by simulations of AGN 
jets powered by the Blandford-Znajek mechanism (see M cKinney] 
( |2006| >) which show an accelerating parabolic base which becomes 
conical near to the point where the plasma reaches equipartition. 

The detailed mechanism behind the production of a relativis- 
tic plasma jet from a black hole is not yet understood fully. It is 
thought that the energy source of the jet comes from the rotation 
of the central black hole and that magnetic field from the accretion 
disc is twisted by the rotation of the black hole creating a Poynting 
flux ( Blandford & Znajek ( 19771). We expect from simulations and 
theory that jets start as Poynting dominated (magnetic field domi- 
nated) and that the particles in the jet are accelerated to a large bulk 
Lorentz factor by the magnetic field until the jet plasma approaches 
equipartition where it reaches a terminal bulk Lorentz factor. The 
accelerating region is thought to be parabolic as indicated by obser- 
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vations fAsada & Nakamura 1 2011), simul ations l |McKinney| ( |2006) ) 
and analytic work (Beskin & Nokhrina (20061), whilst the outer 
regions of the jet are consistent with being approximately conical 
Asada & Nakamura|pOTT),|Hada et al.| ( |2011) , |Sokolovsky et al.| 



201 1 1 and Krichbaum et al. ( 20064) with a slowly decelerating bulk 



Lorentz factor (Mullin & Hardcastle ( 2009 ) and De Young (2010)). 

In this paper we develop our conical jet model from |Pot-| 
ter & Cotter (2012) (hereafter Paper I) to include an accelerating 
parabolic base which transitions to a slowly decelerating conical jet 
after W 5 Rs, motivated by the recent radio observations of M87 and 
consistent with simulations and theory. Our model conserves rela- 
tivistic energy-momentum and particle number throughout the jet 
and includes energy losses to the particles (electrons and positrons, 
hereafter referred to as electrons) via radiative and adiabatic losses. 
We integrate the observed line of sight synchrotron opacity through 
the jet by transforming into the plasma rest frame. We also include 
the effect of inverse-Compton scattering external photons from the 
CMB, an accretion disc, starlight, accretion disc photons scattered 
by the broad line region (BLR) and accretion disc photons repro- 
cessed by a dusty torus and narrow line region (NLR). We calcu- 
late the inverse-Compton emission from these external sources by 
Lorentz transforming into the rest frame of the plasma. This is the 
first time that an inhomogeneous jet model has included an accel- 
erating parabolic base transitioning to a decelerating conical jet. In 
this paper we wish to develop a realistic, physically rigorous jet 
model motivated by and consistent with observations and simula- 
tions. 

A successful jet model should be able to reproduce the simul- 
taneous multi-wavelength observations of blazars across all wave- 
lengths including radio observations, which are usually neglected. 
To test whether our model is capable of producing realistic spec- 
tra we fit it to the simultaneous multi-wavelength spectrum of the 
Compton-dominant FSRQ PKS0227-369 from |Abdo eraLlpOlO) . 

Currently the most popular model of blazars are time- 
dependent fixed radius one-zone models e.g . |Kirk et al.| ( 1998 ), 
|B6ttcher & Chiang| ( |2002| >, |Li & Kusunose| (2000) and |Tsang &| 
|Kirk| ( |2007| i. These models consider the injection of an arbitrary 
electron population into a spherical blob with a constant bulk 
Lorentz factor and calculate the time-dependent emission from syn- 
chrotron, SSC and EC mechanisms. One-zone models best describe 
flaring events observed in the jets, however, they struggle to repro- 
duce radio observations. Investigations into the effects of decelera- 
tion in one-zone models include those by Dermer & Chiang ( |1998| > 
and Bottcher & Principe (2009) where plasmoids are decelerated 
by collisions which accelerate particles. It is important for a realis- 
tic model of blazars to be compatible with radio observations of the 
morphology, core-shift and flat radio spectrum, which are incom- 
patible with one-zone models. 

Investigations of extended inhomogeneous jets include those 



conducted by JBlandford & Konigl |( Ti?7^ , parscTKr1 ( (T9^0l l, |K6nigl 
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Reynolds ( 1982), Mufson et al. ( 1984) and Ghisellini et al. 
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1985 ), Marscher & Gear ( 1 985 ), Kaiser et al. ( 2000 ), Markoff et al 
20001 ), [Sokoiov ^ al.| ( [2004"l >, |Graff et"aTH2008) ), [Pe'er & CaselTa| 
2009| >, |Jamil et al.| ( (20T0| ) and |Vila et aL] ( |201 1| >. Early investiga- 



tions of conical jets were motivated by the original work by |Bland^| 
|ford & K6nigl| ( |1979| l which showed that a conical jet which con- 
serves magnetic energy produces a flat radio spectrum. 

There have been a number of investigations into accelerat- 
ing or decelerating jets including those by Ghisellini & Maraschi 
(l989"l>, |Geo rganopoulos & Marscher] J1998|>, |Spada et al.| ( |2001 

"~ Up 



Georganop oulos & Kazanasj j2003| and Boutelier et al.| q2008f . 
These investigations have included various models for acceleration 



and deceleration in parabolic, conical and cylindrical jets, however, 
no previous work has modelled both the accelerating parabolic base 
and decelerating conical sections of the jet or calculated the line 
of sight synchrotron opacity through the jet. Most investigations 
have variously treated some aspects of the jet physics and emission 
thoroughly whilst neglecting others. The early, pioneering work by 
Ghisellini & Maraschi ( 1989) is the only previous investigation to 
include an accelerating parabolic base transitioning to a conical jet. 
However, this investigation did not consider electron energy losses, 
inverse-Compton scattering of external photons, energy-momentum 
conservation along the jet or a decelerating conical section. In light 
of recent observations and simulations we feel it is valuable to de- 
velop a physically motivated model of jets which includes a thor- 
ough treatment of the relevant physics and which we can use to fit 
to blazar spectra in order to try to answer some of the fundamental 
questions about the blazar population. 

In this paper we will first explain the basic properties of the 
jet including its geometry and variable bulk Lorentz factor. We then 
show how energy-momentum and particle number are conserved 
along the jet. We calculate the observed synchrotron emission and 
inverse-Compton emission from external photons of electrons in the 
jet. We then show the results of fitting the model to the Compton- 
dominant FSRQ type blazar PKS0227. Finally, we discuss the phys- 
ical properties and characteristics of our model fit. 



2 JET MODEL 

We develop the conical jet model presented in Paper I. In this pa- 
per we wish to extend the model to include a variable geometry, a 
variable bulk Lorentz factor, inverse-Compton scattering of external 
photons and adiabatic energy losses. 



2.1 Geometry 

We choose our model to have a parabolic base which transitions to 
a conical jet at a lab frame distance a; con from the black hole. The 
radius of the jet R(x) is given by 



R(x) = R + A pa .rX 



for X < Xcon, 



(1) 



R(x) — -HoH -J 4-par£co^i Xcon) tan(0 O pcning) for X > %con 

(2) 



where x is the lab frame jet length and o 



, is the lab frame 



half opening angle of the jet. We fit the geometry of our model to the 
observed geometry of M87 (Asada & Nakamura (2011)), though 
we leave the opening angle of the conical section of the jet as a free 
parameter since this is likely to vary with the bulk Lorentz factor 
of the jet as ~ — . We assume the geometry of our jet scales 
linearly with black hole mass. This means that the entire geometry 
of our jet is described by only two free parameters (the others are 
set by the observations of M87) the black hole mass M and the lab 
frame conical half-opening angle Opening. 

We do not expect the geometry of the inner region of all AGN 
jets to be the same as that of M87, however, it is very interesting to 
see whether a jet with this geometry is compatible with the obser- 
vations of a blazar. Our model is flexible and so can incorporate any 
length-radius relation if new information on the geometry of other 
AGN jets becomes available in the future. 
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Figure 1, This figure shows the geometry of our jet as used in the model fit 
to PKS0227 where the dimensions are given in terms of the Schwarzschild 
radius of the black hole. The jet starts parabolic with the geometry fixed 
to match observations of M87, the jet transitions to a conical geometry at 
i con = 10 5 -f?s with lab frame half-opening angle ^opening- 



2.2 Acceleration and deceleration 

We set the bulk Lorentz factor of the jet plasma to be a function 
of x. We use the relationship between the bulk Lorentz factor and 
the jet radius predicted by theory (Beskin & Nokhrina (2006)) and 
consistent with simulations (Mc Kinney (2006)) 

TbuikOr) oc x 1/2 . (3) 

Setting the bulk Lorentz factor at the base of the jet to be 
70 and the maximum Lorentz factor at the end of the accelerating 
parabolic region to be 7 max , for x < x con the bulk Lorentz factor 
is given by 



2.3 Energy and particle conservation for the plasma 

In our model the jet plasma starts as magnetically dominated at the 
base, this magnetic energy is converted to bulk kinetic energy of the 
plasma as the jet is accelerated in the parabolic region by a Poynt- 
ing flux ( |McKinney1 ( |2006| l and |Komissarov et al.| < |2009fr ). We do 
not attempt to model the physical mechanism behind the plasma ac- 
celeration since in this investigation we are interested primarily in 
the observed emission of the jet which depends only on the physical 
properties of the plasma. 

We treat the plasma as being an isotropic perfect fluid in its 
instantaneous rest frame along the whole jet. This corresponds to 
the magnetic field being randomly oriented on small scales and 
isotropic and the particles being relativistic with an isotropic ve- 
locity distribution in the rest frame of the plasma. This is an ide- 
alised situation since it requires the magnetic field and electrons to 
isotropise after a section of plasma has been accelerated or decel- 
erated, whilst this is the highest entropy state for the plasma we 
assume the timescale over which the plasma isotropises is short 
compared to the acceleration timescale and this may not be valid. 
However, this assumption of a plasma in which both the magnetic 
field and particle distribution are isotropic and homogeneous in its 
rest frame is used in almost all models of jet emission. 

The magnetic field close to the base of the jet where the plasma 
is Poynting dominated is expected to be approximately parabolic 
( [McKi nney (2006)) and so the electromagnetic stress-energy tensor 
is not that of a perfect fluid, however, in a model of jet emission we 
are only interested in the regions closer to equipartition where high 
energy particles exist or are accelerated. In these regions close to 
or downstream of shocks we expect the magnetic field to be much 
more turbulent and so the approximation of an isotropic tangled 
field is appropriate for these jet regions which contribute signifi- 
cantly to the overall emission. 

We treat the field as tangled on small scales and particles as 
relativistic so that P = |, the energy-momentum tensor in the rest 
frame T'^ u is given by 



7buik(:r) = 7o + 



7max — 70 
^ con 



1/2 

7* ' 



for x < x c 



(4) 



Once the jet has reached equipartition it can no longer signif- 
icantly accelerate the plasma by converting magnetic energy into 
bulk kinetic energy and so the jet stops accelerating and becomes 
conical. The jet now sl owly decelerates jLaing & Bridle ( 2002), 
|Hardcastle et"aT1p005) , pullin & Hardcastle| ( |2009) and |De Young] 
1 2010)) due to interaction with its environment and Compton-drag. 



In this model we will assume the bulk Lorentz factor decreases with 
log(x) . The conical section of the jet starts with a bulk Lorentz fac- 
tor 7 max and decelerates by the end of the jet at x — L to a Lorentz 
factor of 7mi n . Under these assumptions the bulk Lorentz factor is 
given by 



7bulk(a;) = 7max 



7max 7 m i n 



log 



log 



X 

con 



X > X c 



(5) 

The dependence of the deceleration of a jet with jet length is 
not well known so we have assumed this logarithmic decrease. Our 
model is flexible so we can easily change the dependence of the 
bulk Lorentz factor to observations of deceleration if they become 
available. 
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where p is the total energy density in the plasma rest frame and 
p = p e + pB where p e is the particle energy density and ps is the 
energy density of the magnetic field. For a section of plasma which 
is travelling with a bulk Lorentz factor 7buik(£) in the lab frame its 
energy-momentum tensor in the lab frame T^" is given by 



( |7buik(a0 2 p 
37buik(a;) 2 p 


V o 



|7buik(a;) 2 p 
37buik(a;) 2 p 











£ 
3/ 

(7) 

where A M „ is the Lorentz transformation tensor and we have 
assumed that the jet is always relativistic so that j3 w 1. In order 
to satisfy the conservation of energy and momentum along the jet 
our energy-momentum tensor must satisfy V M T M " = 0. In the lab 
frame our jet does not depend on time and quantities only have an 
x dependence so all other derivatives vanish. To conserve the total 
energy-momentum in a slice of width one metre as it propagates 
along the jet we also need to take into account the change in cross- 
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sectional area of the jet (which is invariant under a Lorentz boost in 
the as-direction) as it expands. 



along the jet. The modified electron spectrum contained in one me- 
tre of the jet is given by 



^(■kR 2 (x)T 1v {x)) = J- ^ hu& (x)\R 2 (x)p(x)^ = 0. 

(8) 

We also need to conserve the number of existing energetic 
electrons in the jet as the bulk Lorentz factor of the jet changes, 
so our electron energy density must satisfy the continuity equation. 



AT" (as + das, E' e ) = Norm x N'(x, E' e ). 



(13) 



We now need to ensure conservation of the total energy- 
momentum along the jet. In the accelerating part of the jet we de- 
termine the magnetic field in the section located at x + dx using the 
finite difference formulation of EquationJS] 



V m^ M = £-(*R 2 (x) Pc (x)U(x) 1 ) = 0, (9) 
ox 

where p c is the electron energy density in the rest frame and 
U (as) M is the lab frame four velocity of the plasma. These two con- 
servation equations can be used to calculate the rest frame magnetic 
energy density and electron energy density as the jet accelerates and 
magnetic energy is converted into bulk kinetic energy. We use the 
finite difference formulations of Equations [8] and [9] to evolve our 
magnetic field and electron population along the jet. We do this 
by calculating the synchrotron and inverse-Compton emission of 
the initial electron population, N(x, E e ), in a section of lab frame 
width das. Our code calculates an appropriate width da: adaptively 
such that the shortest radiative lifetime of the population is resolved 
and the radius of the jet does not change by more than 5% over the 
width of the section. We then calculate the resulting electron spec- 
trum as a result of radiative and adiabatic losses over the section 
using Equations [58] and [59] In our code we also allow for the in- 
jection of electrons (with the same initial electron spectrum given 
by Equation|25[l produced from shock acceleration in the section, so 
our modified electron spectrum including energy losses and particle 
injection becomes N'(x, E' e ). 



N'(x,E' e )=N(x,E e ) 



P to t(x,dx',E e ) x Is 
cE e 



+Ai n j(x,dx)E e a e 



,(10) 



where E' e is defined in Equation 57 as the energy E e after 



adiabatic losses have been taken into account and P to t{x, dx' , E e ) 
is the total emitted power from the section of rest frame width dx' 
from Equation [59] The rest frame magnetic energy injected into 
electrons in a slab of rest frame width lm is ttR(x) 2 Ap- m j(x), so 
Ai n j (as, dx) is determined via 



irR 2 (x)A Pill j(x) = / A iuj (x,dx)E, 



1-ct -E./E n 

e ' 



*dE e 



(11) 

In this investigation we will not inject electrons arbitrarily 
along the jet, so we take Ainj (as, das) = 0. We wish to include 
the possibility of an arbitrary injection of electrons in these cal- 
culations, however, since this will allow us to investigate the effect 
of flaring sections of the jet in future work. We use the continuity 
Equation [5] to calculate the overall normalisation of the modified 
electron spectrum in the rest frame of the next section starting at 
x + dx 



■ybuik{x) 2 (p B (x) + p' e (x))%R 2 {x) = 
7buik(z + dx) 2 (pB(x + dx) + p"(x + dx))nR 2 (x + dx), 

(14) 



where the energy densities are defined as 



Pe( X ) = 



f£ Bin E' e N\x,E f e )dE' e 



ttR 2 (x) 



Pb(x) = 



B 2 {x) ^ 
2fi 



p"{x + dx) 



f™ min E' e N"(x + dx,E' e )dE' e 
ttR 2 (x + dx) 



(15) 



In the accelerating part of the jet the initial electron distribution 
at the beginning of the section starting at x + dx is given by 



N(x + dx,E e ) = N"(x + dx,E e ). 



(16) 



We assume that the plasma in the conical part of the jet is al- 
ways in equipartition and the jet is ballistic so no longer suffers 
adiabatic losses as in Paper I. In this case the bulk kinetic energy 
of the jet is converted into the acceleration of electrons at shocks 
and amplification of the magnetic field as the jet decelerates and so 
we use Equation[8]to determine the amount of energy injected into 
electrons through shock acceleration in a slab of rest frame width 
lm Ap dcc (as + das)-7r i?(as + dx) 2 . 

7buik(x) 2 (pb(x) + p' c (x))nR 2 (x) = -y hulk {x + dx) 2 ... 
x (/Ob (a; + das) + p"(x + dx) + Ap dc c{x + dx))nR 2 (x + dx), 

(17) 

where we ensure the jet remains in equipartition in the conical 
section via the constraint 

Pb(x + dx) = p"(x + dx) + Ap dcc (as + das). (18) 

The electron spectrum at as + das in the decelerating part of the 
jet is given by 



N(x + dx,E e ) = N"(x + dx,E e ) + C(x + dx)E~ a e 



a -E e /E n 



Apdec(as+das)7ri?(as+das) = C(x+dx) / E e a e 

' Em 



(19) 

.., VK,.n* d/ ,; 



(20) 



Norm ; 



7buik(as)-R 2 (as) 



7buik(a- + dx)R?(x + dx) ' 



(12) 



This takes into account the effect of differential length con- 
traction on the number density of electrons as the population moves 



3 PHYSICAL CONDITIONS AT THE BASE OF THE JET 

We wish to calculate the magnetic field strength and electron pop- 
ulation in the rest frame in terms of the physical parameters of the 
jet. For a relativistic jet with total kinetic luminosity of Wj, — is 
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equal to the total energy in a section of plasma of width one metre 
in the x direction in the lab frame Et . 



Ej = T m (x = 0)dV = Utp(x = 0)nR 2 , 



p(x = 0) = 



3Ej 



(21) 



(22) 



p(x = 0) = p B (x = 0) + p e (x = 0) , A cqu i (a;) 



Pb{x) 

Pe(x) 



p B (x = 0) = 



-PQ 

2/xq ' 



(23) 
(24) 



We assume that the high energy electrons have been acceler- 
ated through sho cks ptell| ( |1978 1) or magnetic reconnection ( |Zeni-| 
|tani & Hoshino| ( |2001^ "and |Zenitani & Hoshinolj2007> ). Our code 
is flexible and allows us to inject an arbitrary electron population 
at different distances along the jet so we are not tied to a specific 
mechanism for particle acceleration. In this paper we characterise 
the spectrum as being a power law with index a with an exponential 
cutoff at energy E max and with a minimum energy U m in. A power 
law electron distribution with an exponential cutoff is an approxi- 
mation to the steady state solution to the diffusion loss equation for 
a shock including radiative energy losses (see for example Bregman 
( |1985| >). The number of electrons contained in a slab of width one 
metre in the rest frame of the plasma is given by 



N e (x = 0, E e ) = A{x = 0)E~ 



(25) 



In general our electron distribution will evolve due to energy 
losses and additional injection of electrons to take a more general 
form where A(x, E e ) becomes a function of electron energy and 
length along the jet. Our initial value of electron energy density is 
therefore 



^Dopplcr 



Or)' 



^Dopplcr 



(x) = 



l 



TbuUt 0*0(1 



P(x) COS(0 o bs)) ' 

(27) 



cos(6> 



fi{x) 



P(x) cos(6*o 



7buik(a;) 



(28) 

We track photons emitted from interior sections of the jet as 
they move through outer sections by transforming to the rest frame 
of each section. For a given observed lab frame frequency v all 
the contributing synchrotron photons from different sections with 
different bulk Lorentz factors will have the same frequency v' and 
emission angle 9' obs in the rest frame of the section at x. This means 
that to calculate the emission at an observed frequency v we need 
to calculate the synchrotron emission at different frequencies v' in 
sections of plasma with different bulk Lorentz factors. To calculate 
the line of sight opacity to a section we need to integrate the opacity 
at rest frame frequency v' across the rest frame path lengths of all 
the outer sections. 

We use the result from Paper I for the distance travelled 
through the plasma by a synchrotron photon in the rest frame 



dr' = 7buikda; 



COS(# 



■/*(*: 



(29) 



where dr' is the total rest frame distance travelled by the pho- 
ton. We use the equations for the synchrotron emission and opacity 
k v (y, x) from Paper I with the exception of Equation 23 in Paper I 
which contains an erroneous extra factor of 2 in the numerator. Us- 
ing the formulae for the synchrotron opacity from Paper I we find 
the observed optical depth to a section to be given by 



Ttot{v,x)= I k v {v'(x), x)dr'. (30) 

J X 

The observed synchrotron emission from the jet is given by 



Pe(x = 0) = 



E N.( X = 0,E e ) 
■kR 2 



(26) 



These equations specify the evolution of the plasma along 
the jet conserving energy-momentum and electron number density. 
We now turn our attention to calculating the observed synchrotron 
emission from the jet. 



vF v = ^fooppicr^) 4 ^ P v (x,dx,v')e Ttot{ - v -^ 



(31) 



where the length of an emitting section in the rest frame of 
the plasma dx is related by a Lorentz contraction to the lab frame 
section length dx. 



4 SYNCHROTRON EMISSION 

In this section we shall generalise the calculation of synchrotron 
emission from Paper I to the case of a jet with a variable bulk 
Lorentz factor. In the lab frame we wish to calculate the observed 
luminosity at a given observation angle of the jet. To do this we 
will calculate the emission and absorption of synchrotron radiation 
corresponding to an observed lab frame frequency and observation 
angle as it propagates along the jet by Lorentz transforming into 
the rest frame at every section of the jet. Observed synchrotron ra- 
diation at a frequency v and angle 8 obB corresponds to an emitted 
frequency v' and angle 6' ohs in the rest frame of a section of plasma 
with bulk Lorentz factor 7buik given by a Lorentz transformation 



5 INVERSE-COMPTON EMISSION 
5.1 SSC emission 

To calculate the SSC emission we use the result of the emission 
produced in the plasma rest frame from Paper I. Since the high en- 
ergy emission is optically thin the only modification to the earlier 
result is that we need to Doppler-boost the emission from each sec- 
tion individually. We use Equation 38 from Paper I to calculate the 
inverse-Compton power emitted by a section in its rest frame per 
log bin of rest frame frequency (corresponding to v' Fl). 

vF v = ^ v'(x)F' v {v'{x),x, dx)(5 Do ppi cr Qr) 4 . (32) 
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5.2 Inverse-Compton scattering of external photons 

In this work we consider the six most likely sources of external 
photons to the jet: the CMB, radiation directly from the accretion 
disc, starlight and radiation from the accretion disc scattered by free 
electrons in the BLR and reprocessed by a dusty torus and NLR. 
Inverse-Compton scattering of external photons has been investi- 
gated by a number of auth ors including |Dermer et al.H1992),|Der- 



|mer & Schlickeiser|( [T993l , |Sikora et al.| ( |1994| >, |B6ttcher & Bloom 
pOOOl i and Hardcast e & Croston| ( |20TT) . 



The full calculation for the inverse-Compton emission of an 
axisymmetric external photon field by an isotropic electron distribu- 
tion takes the form of a seven dimensional integral. In order to make 
the problem computationally tractable we make a number of sim- 
plifying approximations. Firstly we assume that in a collision with 
an electron the scattered photon travels with a velocity parallel to 
the electron's initial velocity. This approximation is commonly used 
when considering inverse-Compton scattering of photons ( Dermer 
|et al.| p009 l) and is justified because the scattered photon's energy 
is a steep function of its scattering angle with significant scatter- 
ing energies occurring for photons travelling close to the electron's 
initial velocity, due to conservation of momentum. 

For a jet with a significant bulk Lorentz factor radiation which 
is isotropic relative to the black hole is Doppler-boosted in the rest 
frame of the jet plasma, such as the CMB, starlight, BLR, NLR 
and dusty torus photons. The radiation is Doppler-boosted in the 
plasma rest frame so that it is beamed along the jet axis in the 
opposite direction to the jet material with a characteristic opening 
angle ~ — - . This means that for large values of 7buik almost 
all external photons travel close to the jet axis. We shall make the 
approximation that all the external photons travel in the negative 
^-direction to reduce the dimensionality of the integral to 3. This 
approximation is strictly valid only in the limit of large 7buik, how- 
ever, the scattered energy of the photons depends on the incident 
photon direction through oc (1 + /3 cos(6>)) so for a small change 
in opening angle — - — the scattered photon energy will not change 
significantly. The total number density of external photons is there- 
fore increased by a single Doppler factor and their energy is also 
increased by a Doppler factor in the rest frame of the plasma (since 
we have effectively integrated over solid angle). 

In the case of external photons travelling directly from the ac- 
cretion disc we make similar approximations. We approximate the 
scattered photon velocity to be parallel to the initial electron ve- 
locity. We also approximate the external accretion disc photons to 
travel in the positive a;-direction, this is justified because the bright- 
est part of the accretion disc is closest to the black hole (at a few 
Schwarzschild radii) and so the base of the jet is not illuminated 
isotropically by the accretion disc. The most energetic photons illu- 
minate the disc from behind in the frame of the black hole. 

Unlike in the case of photons from the CMB, starlight, BLR, 
NLR and dusty torus, photons directly from the accretion disc are 
de-boosted, so Lorentz transforming into the plasma rest frame re- 
sults in the photon angles diverging away from a>axis. This means 
that our approximation is only justified in the case that the radius 
of the jet is smaller than the distance from the accretion disc. If 
this is not the case then the approximation results in an underesti- 
mate of the inverse-Compton emission from accretion disc photons, 
however, neither the Doppler factor or scattering cross section de- 
pend very steeply on the angle of the photons from the a>axis for 
the case of photons illuminating the jet from behind so the results 
should not differ significantly from a more detailed treatment (see 
for example Der mer et al.|p009| ). Since most previous investiga- 



tions have assumed de-boosted accretion disc photons are not rel- 
evant to inverse-Compton emission it is interesting to see whether 
this emission is significant. 

With these approximations we can use a modified version of 
the inverse-Compton scattering formula from paper I (Equation 36) 
to calculate the inverse-Compton emission from external photon 
sources. Using these approximations we have eliminated the inte- 
gration over 8. Instead we only consider scattering from electrons 
travelling at an angle 8 where this angle corresponds to photons 
scattered in the lab frame at the observation angle of the jet # obs . 
This angle 9 is easily related to the observation angle of the jet via a 
Lorentz transformation which in the case of CMB, starlight, BLR, 
NLR or dusty torus photons we calculate to be. 



(33) 



cos(6> obB ) - /3(x) 
1 - fi{x) cos(6» obs ) J ' 
and for the case of external photons directly from the accretion 



disc 



cos(6> obs ) - /3(x) 
1 - P(x) cos(6» obs ) 



(34) 



Fixing the value of 8 using one of Equations[33]or[34]we find 
the inverse-Compton emission of external photons to be 



weight (-E e , .E 7 , 0, <j>2,x) = N e (E e ,x) 



xn 1 (E 1 ,8, x).c(l — P e cos 8)ir \ sin 02 1- 



(35) 



Power(W0 = E -<s 

catt .weight(_/? e , -/?-y, 8, (j>2-, x} 

xd£ e d£ 7 d0 2 ds'. (36) 

where /3 e is the electron velocity divided by c. We now proceed 
to calculate the number density of external photons in the plasma 
rest frame for the six sources. In the following subsections quanti- 
ties defined in the lab frame of the jet are primed and quantities in 
the rest frame of the jet plasma are not primed. 

5.3 CMB 

We use the blackbody distribution to find ra 7 for the CMB at redshift 
2 withT = 2.73(1 + z)K. 



n'Jv) 



47T 



2v'- 



(37) 



We average the Doppler factor over solid angle to find the 
mean Doppler factor for a photon emitted isotropically in the black 
hole rest frame which is 7 bu ik- The frequency of the photons and 
their arrival rate are both boosted by a Doppler factor. The number 
density per unit frequency of photons is therefore unaltered since 
the increase in arrival rate is compensated by the stretching of the 
frequency bin occupied by the photons. 

6i> pp(x)n\(i/')du' = rij(u)du, Su opP (x) = 7buik, (38) 

n 7 (» = 5 D op P (3;)n 7 ( - — ^j-t) j % , , , n^iy) = n 7 (i/ 

\ovopp\x) ) 0v O pp(x)dv' 

(39) 
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From these equations we find that the total energy density of 
CMB photons in the plasma rest frame (7 p hot = / hvn~f{y)Av 
is proportional to bulk Lorentz factor squared as we expect from 
Equation|7] 



5.4 Thin Accretion Disc 



We consider a thin accretion disc ( Shakura & Sun yaev| ( |1973) ) 
around a black hole of mass M accreting at a rate M. Since we 
do not know the likely black hole spin parameters of these AGN we 
will conservatively assume a slowly rotating black hole. We treat 
the thin disc as a radiating blackbody. 



3GMM 4 
dL = — — — — dr = 2aT 2nrdr, 
Ir 2 - 



T = 



( 3GMM\ 



L v = 



\ 8anr 3 J 
4whri/' 3 

dr, n 



6GM 



(40) 



(41) 



(42) 



We make the approximation that the accretion disc luminosity 
is emitted isotropically and consider the photon density at a radius 
R from the disc in the lab frame. 



4nR 2 hv'c 



(43) 



The photon number density and frequency both decrease by a 
Doppler factor when they are Lorentz transformed into the plasma 
rest frame. Again we use the value of the Doppler factor averaged 
over solid angle 



n^(v) = n' 7 (i/ST> ap p(x)), foop P (:r) = 7 bu ik- 



/ \ Pl[t = 40kpc)(40kpc) _ 3 
Pi(r) = £ — ^ j-^ -(Wpc ) for r>40kpc, 

(46) 

where the luminosity density is assumed to be constant within 
the central lpc of the galaxy. We calculate the total enclosed stel- 
lar mass in the core of the galaxy corresponding to this luminosity 
density to be 4.0 x 10 12 Mq using a mass to light ratio of 10, this 
is consistent with the expected stellar mass in the core of a large 
elliptical galaxy. We calculate the energy density of starlight ob- 
served by the jet plasma at a distance x along the jet (/9 s tar(a0) by 
the integral 



lOOkpc 



dr, (47) 



(44) 



where we have started the integral at the position of the jet sec- 
tion, x, since the starlight from regions r < x will illuminate the jet 
from behind and so will be strongly Doppler-deboosted. We have 
also made use of the spherical symmetry of the problem to sim- 
plify the integral, since all observers inside a spherically symmet- 
ric emitting surface see the same isotropic luminosity and surface 
brightness (due to Newton's shell theorem) we are free to calculate 
the observed energy density of an observer at the origin to simplify 
the integral since this is the same as the luminosity observed by the 
plasma at r = x. Using this integral we find the energy density of 
starlight at a distance of lpc from the centre of the galaxy to be 
~ 2 x 10~ 12 Jm~ 3 . This is in agreement with the value quoted in 
|Hardcastle & Croston| ( |20iT) for the energy density in starlight at 
the centre of an elliptical galaxy. 

We assume that the starlight is predominantly due to K-type 
stars which typically dominate the spectra of elliptical galaxies. We 
choose the photon distribution to be that of a 5000K blackbody 
since this corresponds to the surface temperature of a typical K-type 
star. In practice it is difficult to accurately model the interstellar ra- 
diation field at the centre of a powerful FSRQ since the detailed 
properties of the stellar population will vary substantially between 
galaxies. We include starlight in our model for completeness and to 
investigate whether there is evidence of inverse-Compton scattering 
of starlight from fitting to the spectrum of PKS0227. 



5.5 Starlight 

Inverse-Compton scattering of starlight by jets has been investi- 
gated previously by Hardcastle & Croston (2011 1. In this work we 
are fitting to a powerful Compton-dominant blazar so we expect its 
host galaxy to be a giant elliptical. Work by Lauer et al. ( 1995} and 
|Faber et HT] J1997) found that the luminosity density at the centre 
of giant ellipticals is considerably lower than that of smaller ellipti- 
cal or disk galaxies with the luminosity density at 30pc in the most 
massive ellipticals found to be typically between 2 — lOZ/0pc~ 3 . 
These large ellipticals have relatively flat cores over which the lu- 
minosity density goes as pL oc r a with —1.3 < a < — 1.0 ( jLauerj 
|et al.| ( [1995b), at larger distances th e luminosity density decreases as 
pL oc r~ |Tremaine et al. 1 1994 1). In M87 the flat core with a lu- 
minosity density power law index a = —1.3 extends out to 40kpc 
( |Cohen & R yzhov ( 1997 1). Using the typical values above to model 
the seed photon population due to starlight along the jet we choose 
a luminosity density given by 

Pl{t) = 5£ e(y ,2 ( Wpc -3) for i pc > r <40kpc, 

(45) 



5.6 Accretion disc photons scattered by the broad line region 

To calculate the number density of accretion disc photons scattered 
by the broad line region in the lab frame we use the result from 
IDermer et al~1 l |2"009l 



n'JE'\ 



t (R)a T L ! , cc (E^ 
A-kRcE' 



(48) 



where n® LB (R) is the free electron number density. The au- 
thors calculated this result assuming that accretion disc photons are 
Thomson scattered by free electrons in the BLR. The result assumes 
that the scattered radiation field is isotropic in the lab frame so the 
photon number density and energy in the plasma rest frame are re- 
lated to that in the lab frame by a Lorentz transformation. 



EL 



<5Dopp(2) 



8uop P {x) = 7buii 



(49) 



where we have used the value of the Doppler factor averaged 
over solid angle. We assum e the BLR extends out to .Rblr = 



( 2005 )), which we take 
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to be 0.5pc for the most powerful FSRQ type blazars such as 
PKS0227. To calculate the value of LR we assume a constant 
number density of electrons out to -Rblr and a covering fraction for 
the BLR, 77blr =0.1 ( |Ghisellini & Tavecc hio (2008) and |Sikora| 
|et al.| ( [2"0"09] >). The covering fraction is given approximately by 



T]BLR 



BLR 

-a n9 T 47Tfl 2 dfl » n^ LR a Tj RBLR, (50) 



For our chosen covering fraction and BLR radius this gives 
n BLR = 1.0 x lO^m -3 . 



5.7 Dusty Torus and Narrow Line Region 

We now calculate the number density of photons scattered by 
clumps of dust in the torus and narrow line region (NLR) surround- 
ing the black hole. At a distance x from the black hole the surface 
temperature of the dust facing towards the accretion disc, Td us t, is 
given by ( |Nenkova et al.H2008) ) 



1500 



K 



T dust = 640 



x d J 



x < 9x d , 



x > 9xd- 



(51) 



where x d = 0.4(Z, acc /10 37 VT) 1/2 pc is the dust sublimation 
distance and L acc is the total accretion disc luminosity. We as- 
sume a clumpy dust cloud distribution where the number density 
of clouds is given by n c i ou d x ~ q with 1 < q < 2 follow- 
ing Nenkova et al. (2008). We also assume that the distribution of 
cloud sizes is independent of x. The hot torus absorbs and repro- 
cesses a fraction (the covering fraction) t)t of the incident accre- 
tion disc radiation. We use a covering fraction of rjr = 0.1 taken 
from the ratio of accretion disc to dusty torus emission in Nenkova 
|eTaIJ ( [2008l l. We further assume that this hot torus extends from the 
dust sublimation distance Xd out to ~ 10xd as indicated by obser- 
vations ( |Elitzur| ( [2006| > and |Nenkova et al.| ( [2008l >). We assume that 
the NLR has a covering fraction ??nlr = 0.07 ( |Mor et alT] ( |2009) ) 
and extends from the outer radius of the torus (x = lQxd) out to a 
distance :r outer , typically around lkpc (Be nnert et al.|(200"6| l). 

We now calculate the number density of photons from this 
dust distribution. Taking an average cross-sectional area of a cloud 
(Jdoud independent of x, the cloud number density and covering 
fraction we can calculate the emissivity of the dust as a function 
of distance from the black hole. The average emissivity of the dust 
j(x) at a distance x will be proportional to the number density of 
clouds at x, so j(x) = jox~ q 



use Newton's shell theorem to calculate the energy density of radi- 
ation inside an emitting shell of dust. We assume clouds re-emit the 
absorbed accretion disc radiation as a black body with a tempera- 
ture given by Equation|5T](this is the correct dust temperature since 
we are only interested in the emission from the surface of clouds 
facing back towards the jet and black hole). In this case the number 
density of photons from dust outside a distance x in the lab frame 
is given by 



(53) 



where j(x, v') is the appropriate blackbody distribution nor- 
malised to the correct emissivity. 

., _ TTj X-" 2hv ri 1 

J[X,V ' ~ oT d ™t(x) 4 ' c 2 ■ e ^'/fcBT du8t W_i- (M) 

The number density of photons in the plasma rest frame is then 
given by 



n v {x, v) — n I x, — - — ) , foopp = 7buik(a;). 

V ODopp / 



(55) 



In this work we take the value of q — 1 for the cloud distribu- 
tion as suggested by |Mor e~ al. (2009) from fitting to observations. 

5.8 Observed external Compton emission 

The observed inverse-Compton emission from the external photons 
is related to the emission in the rest frame by 

vF v = ^ v {x)F' v {y' {x) , x, dz)fe ppi cr (:r) 4 . (56) 



6 PARTICLE ACCELERATION AND ENERGY LOSSES 

We take into account radiative and adiabatic losses to the electron 
population as it traverses along the jet. We use Equation[59]to take 
into account the effects of radiative losses to the population due to 
both synchrotron and inverse-Compton emission. In addition to ra- 
diative losses in this investigation we will also consider adiabatic 
losses in the accelerating parabolic section of the jet. We assume 
the conical section of the jet is ballistic and so does not suffer adi- 
abatic losses as in Paper I. We use the formula from Kaiser] < [2006) > 
to include adiabatic losses to a relativistic electron population. 



E' = E 



= E 



(57) 



r\T Ls.cc = const 



[10x d 



n c ioud(2:)o-cioud47rx da: 



jo / 4nx 2 x q dx. (52) 



From this equation we can calculate jo for both the torus and 
NLR (with the appropriate torus and NLR quantities interchanged). 
Since the bulk Lorentz factor of the jet is large we only consider 
emission from spherical shells at a radius greater than the distance 
of the jet section from the black hole, since radiation emitted from 
behind the jet is Doppler-deboosted in the rest frame of the jet. We 



N'(E') = N(E), 



(58) 



where E is the initial energy of an electron, V is the initial vol- 
ume of the section and R the initial radius, primed quantities are 
those after expansion. Considering the electron distribution con- 
tained in a slab of width one metre in the rest frame we use the 
following equation to take into account radiative energy losses to 
the electron population as it travels along the jet. 



N e (E e ,x + 7buikda;') = N e (E e ,x) 



P to t(x,dx',E e ) x Is 



cE e 



(59) 
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where P tot is the total power emitted by electrons of energy 
E e in a section of rest frame jet width dx' due to synchrotron and 
inverse-Compton emission. The energy losses for the synchrotron 
and SSC emission are the same as those calculated in Paper I. 
We have simplified the calculation of the inverse-Compton scat- 
tering of external photons by considering only the scattering events 
which lead to photons travelling at the observation angle. To cal- 
culate the electron energy losses from this anisotropic photon field 
we must take into account the average energy lost from isotropic 
electrons encountering a narrow beam of external radiation. In this 
case we can use the formula in Paper I for the total power emitted 
by an isotropic electron distribution scattering an isotropic photon 
field using the value of n 1 calculated above for the external photon 
fields. After radiative losses have been taken into account we then 
apply Equation [58] to the population to take into account adiabatic 
losses in the accelerating parabolic section. 



7 ACCELERATION OF ELECTRONS IN THE 

DECELERATING CONICAL PART OF THE JET 

Observations of AGN jets show evidence for in situ acceleration oc- 
curring along the jet. High energy optical synchrotron electrons are 
observed at distances beyond that expected from their synchrotron 
lifetime (see Meisenheimer et al. ( 1997} and |Jester et al.| ( [2001| >). 
Observations also indicate that jet material decelerates away from 
the base though the jet remains relativistic until the hot spot l Laingj 
|& Bridle|p002) , |Hardcastle et al.|p005] l and |De Young] ( |2010| >). In 
our model shocks in the jet convert bulk kinetic energy into elec- 
tron energy in the rest frame of the plasma as the jet decelerates 
once it has become conical. This is a natural result of conservation 
of energy-momentum in our model since a decelerating jet loses ki- 
netic energy and so must gain internal energy (or do work on its 
environment) to conserve energy. We have assumed that the excess 
bulk energy is converted into both particle acceleration and ampli- 
fication of the magnetic field to maintain equipartition throughout 
the conical section of the jet. 

In this work we simulate the effect of acceleration of electrons 
by the injection of a fixed power law distribution of electrons con- 
tinuously along the jet which conserves energy-momentum in the 
decelerating part of the jet. Observations of AGN jets show a contin- 
uous variation in optical spectral index along the jet with evidence 
for continuous acceleration taking place along the jet ( Jest er et al.| 
( |2001| >) and so converting bulk kinetic energy to internal electron 
energy through deceleration is a natural way in which to explain 
these observations. The acceleration of these electrons is thought to 
be caused by shock acceleration where the shocks originate from a 
variety of mechanisms including the collision of shells of plasma 
travelling at different speeds as the jet decelerates (|Spada et al.| 
P001| >), shear acceleration ( Rie ger & Mannheim| (2002 1) or ex- 
ternal pressure driven recollimation shocks (Komissa rov & Fallej 
( |1998| >). Variations of models in which in situ acceleration of elec- 
trons occurs are popular and have been investigated by many au- 
thors including [Kaiser et al.lj2000^ , |Spada et al.| < |2001| > and |Jamil| 
|et al.|p0T0| l. This is the first time, however, that conversion of bulk 
kinetic energy into accelerating electrons through deceleration has 
been considered self-consistently in a continuous conical jet with a 
variable Lorentz factor. 

We assume that the injected electrons have the same spectrum 
as the initial electron population (Equation |25[ >. In this investigation 
we will assume that the steady-state shock distribution is the same 
for all shocks along the jet. In practice the power law exponent and 



maximum energy cutoff will depend on the detailed properties of 
the shock and plasma (Bell ( 1978 1) but this is beyond the scope of 
our current investigation. 

We calculate the evolution of the electron distribution along 
the jet by taking into account radiative energy losses on the electron 
population (and adiabatic losses on the remaining electrons in the 
parabolic section). We determine the energy density of electrons 
injected into the plasma due to deceleration by Equation|17| 



8 EQUIPARTITION CONSIDERATIONS 

Observations indicate that on large scales the jet plasma is close 
to equipartition (Croston et al. (2005)). Observations of brightness 
temperatures of the cores of jets also indicate that the jet plasma 
is close to equipartition when sources are in a quiescent state, 
however, in some cases shocks are believed to produce regions 
where the plasma becomes particle dominated associated with flar- 
ing ( [Homan et al.| l |2006} ). It is thought that jets are produced via 
the interaction of electromagnetic fields with a rotating black hole 
(Blandford & Znajek ( 1977)). Current models of jet production start 
with a magnetically dominated region which accelerates particles to 
form the jets we observe. Simulations indicate that the jet becomes 
conical when the plasma nears equipartition ( McKinney ( 2006) and 
|Komis sarov et al. ( 2009 1). This makes physical sense since the mag- 
netic field has done work accelerating the electron plasma and in 
equipartition the remaining magnetic energy is not large enough to 
accelerate the jet substantially. In this model we force the jet plasma 
to be in equipartition after the point at which the jet becomes coni- 
cal by the injection of shock accelerated electrons. 

For a given total jet power we expect the synchrotron and SSC 
emission to be brightest when the jet is close to equipartition from 
a reversal of the minimum energy argument for synchrotron emit- 
ting lobes (Burbidge ( 1959)). If the accelerating region is magneti- 
cally dominated we therefore expect the synchrotron emission to be 
brightest at the point where the jet becomes conical since the bulk 
Lorentz factor, jet power and magnetic field will be largest at this 
point and decrease along the conical section. 



9 CODE 

We use the equations given in the preceding sections and Paper I 
to calculate the emission from the electrons along the jet. In order 
to make the problem computationally tractable we must discretise 
the electron population, sections of the jet and the particle accel- 
eration resulting from the deceleration. Our code uses logarithmic 
bins of electron energy where the bin of energy E e has a width 
dE e — O.SE e , so the total electron population for PKS0227 spans 
approximately 40 energy bins and for high frequency peaked BL 
Lacs this increases to approximately 60 bins. The jet sections are 
determined adaptively so that either the shortest radiative lifetime 
is resolved or the jet radius does not increase by more than 5%, 
this results in jets typically being composed from 200-600 cylindri- 
cal sections. At each decelerating section the equipartition fraction 
of the jet is calculated and additional electrons ( Apdcc in Equation 
|17} are only injected if this equipartition fraction differs from 1 by 
more than 10% so that we are not continually limited to resolving 
the shortest electron radiative lifetime. 

We have chosen the above values for discretisation so that in- 
creasing the resolution/number of bins does not noticeably affect 
the spectrum produced from the jet, whilst this allows our code to 
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produce a spectrum for a jet in 15-30 seconds running on a single 
processor core. 



10 RESULTS 

In this paper we have set out to develop a physically realistic jet 
model which is consistent with observations and simulations and in- 
cludes all the relevant physics to accurately predict its emission. In 
order to test whether the model can potentially describe real blazars 
we shall try to fit to the spectrum of the most Compton-dominant 
FSRQ from flAbdo et al.| pOlO) ) PKS0227-369. We choose a 
Compton-dominant blazar since in Paper I we showed that a sim- 
pler version of this model was able to reproduce the spectrum of 
BL Lacertae, a relatively weak powered object. Compton-dominant 
blazars are thought to represent high-powered blazars (Foss ati et al.| 
( |1998| >) so it is interesting to see if our model is capable of reproduc- 
ing the spectra of the two extremes of the blazar population. We as- 
sume a standard ACDM cosmology ( Komatsu et al. 1 201 to cal- 
culate the distance to PKS0227 which has a redshift of z = 2.115 
( |Abdo et al.| ( |20T0l ». 

We will first show the results of fitting our model to the data 
using a transition region at different distances along the jet at large 
distances > lOpc, within the dusty torus ~pc and within the BLR 
< 0.5pc. One of the most important current questions about jets is 
where the gamma-ray emission originates from. Using our realistic 
continuous jet model we attempt to fit the quiescent jet emission 
of PKS0227 and try to constrain the location of the quiescent jet 
emission. 



10.1 Fit 1 - A transition region at large distances 

We first attempt to fit the spectrum using a transition region at large 
distances outside of the dusty torus and BLR as implied by the ge- 
ometry of M87 where the jet transitions from parabolic to conical 
at 10 5 i? s , corresponding to a distance of 66pc. 

In Figure[2]we show the results of fitting the model by eye to 
the simultaneous multi-wavelength spectrum of PKS0227. We see 
that the model fits very well to the observations across all wave- 
lengths including radio data points with a transition region at a lab 
frame distance of 34pc and a jet radius of 0.64pc. In Table [l] we 
show the model parameters for this fit to PKS0227. We find that the 
fit to the observations requires a high power (Wj = 1.2 x 10 39 WO, 
high bulk Lorentz factor jet (7 max = 45), observed at a small an- 
gle to the line of sight (# bs = 0.8°). This is consistent with our 
expectations of Compton-dominant FSRQs being FRII type AGN 
observed close to the jet axis (Urry & Padovaru] ( [1995) ). We find 
that the jet has a maximum bulk Lorentz factor of 45 and the con- 
ical section of the jet decelerates to a bulk Lorentz factor of 35 by 
x = 2 x 10 21 m. For these parameters the magnetic field strength 
in the plasma rest frame at the transition region is relatively low, 
2.03 x 10 T, required by the low synchrotron peak frequency and 
Compton-dominance of the blazar. We find a value of the electron 
spectral index a = 1.9 close to the value expected from standard 
diffusive shock acceleration. 

At these large distances the Compton-dominance is due to 
inverse-Compton scattering of CMB seed photons with a small con- 
tribution by NLR photons at the highest gamma-ray energies (see 
Figure[3k). In Figure|4]we show the energy density of the different 
radiation fields measured in the plasma rest frame at different dis- 
tances along the jet for this fit. We see that the NLR dominates at 



the transition region but the CMB is the dominant photon source at 
large distances along the jet. 

We include starlight seed photons in our model, however, we 
find that inverse-Compton emission from the CMB dominates the 
emission from starlight, as shown in Figure[3k, except at the highest 
energies. This is because the energy density in the CMB increases as 
oc (1 + z) 4 , so whilst the energy density in starlight is significantly 
larger at z = 0, we find that the CMB photons are more signifi- 
cant for high-redshift Compton-dominant blazars such as PKS0227 
(see Figure ffl. The inverse-Compton scattering of starlight is also 
reduced relative to that of the CMB due to the decreased Klein- 
Nishina cross section for the higher energy starlight photons. 

We find that accretion disc photons directly from the disc are 
the dominant photon energy density only close to the base of the 
jet (within a few Rs of the black hole) and quickly become dom- 
inated by other photon sources due to the Doppler-deboosting and 
the 4j dependence of the photon number density (see Figure|4j. We 
find that for our fit to PKS0227 the Compton-dominance is due to 
inverse-Compton scattering of CMB and NLR photons in the con- 
ical jet region, as shown in Figure [3k. Previous investigations tend 
to ignore the contribution of CMB photons to the inverse-Compton 
emission, this is because in a BL Lac type object with a low bulk 
Lorentz factor jet at low redshift the CMB is dominated by SSC 
photons. We find that due to the high bulk Lorentz factor and high 
redshift of our model fit the CMB becomes the dominant photon 
source at large distances along the jet, > a; ou tor, when taking into 
account the Doppler-boosting of the CMB into the rest frame of the 
plasma and including the redshift dependence of the CMB number 
density. 

This is an interesting result since it has implications for the 
Cosmic evolution of Compton-dominant blazars. If their Compton- 
dominance is due to inverse-Compton scattering of CMB photons 
then we might expect to see a redshift dependence of the Compton- 
dominance of blazars, with increasing Compton-dominance at 
higher redshift (since the energy in the CMB scales as (1 + z) 4 ). 
This could explain why most Compton-dominant blazars are ob- 
served at high redshifts. We might also expect to see rarer high 
powered, high Lorentz factor blazars at large redshift due to the 
increase in comoving volume and due to the larger number of ma- 
jor mergers in elliptical galaxies which host FRII sources. In prac- 
tice it may be difficult to disentangle all these evolutionary effects, 
however, we might expect that high Lorentz factor blazars at low 
redshift are less Compton-dominant than the same sources at high 
redshift where the CMB has a higher blackbody temperature. 

Figure |3J) shows the total emission from different sections of 
the jet. The synchrotron emission is dominated by the transition 
region of the jet due to its high bulk Lorentz factor and higher mag- 
netic field strength than outer parts of the jet. The low frequency 
radio emission originates from farther along the jet where the jet 
radius is larger. The inverse-Compton emission is dominated by the 
transition region and outer parts of the jet where electrons are ac- 
celerated through shocks in the decelerating part of the jet and ra- 
diate this energy slowly, primarily through scattering NLR photons 
at distances x < x out cr and CMB photons at larger distances. 

In Figure [5] we show the evolution of the electron population 
in a slab as it travels along the jet in fit 1. A small amount of energy 
is injected into electrons at the base of the jet which quickly cool 
both radiatively and adiabatically from high to low energies due to 
the large magnetic field strength. The population is replenished by 
the shock at the base of the conical section located at 10 5 R a = 
1.08 x 10 18 m where the plasma is in equipartition. The population 
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Figure 2. This figure shows the results of fitting the model to the SED of PKS0227. The model fits the observations very well across all wavelengths including 
radio observations. We find the model fit requires a high jet power, high Lorentz factor jet observed at a small angle to the jet axis. We find that the inverse- 
Compton peak is due to scattering of Doppler-boosted high-redshift CMB with a small contribution from scattering NLR photons at the highest gamma-ray 
energies. 




eV eV 

(a) The fit to PKS0227 showing the different contributions to the inverse- (b) The fit to PKS0227 showing the emission from different distances 
Compton emission. along the jet 

Figure 3. This figure shows the different inverse-Compton components in Figure a and the jet emission from different distances in Figure b, using a transition 
region at large distances. Figure a shows that the inverse-Compton emission is produced by scattering CMB photons with a small contribution from NLR at the 
highest gamma-ray energies. Figure b shows that the optically thin synchrotron emission is dominated by the transition region, with the radio emission coming 
from larger distances as we expect. The inverse-Compton emission originates at large distances from the central black hole from scattering CMB photons. 



then suffers radiative losses as it travels away from the shock though 
less severely than at the base of the jet. 



10.2 Fit 2 - A transition region within the dusty torus 

We now try to fit the spectrum using a transition region within the 
dusty torus but outside of the BLR. An emission region within the 



dusty torus has been investigated previously using one and two- 
zone models by jiBlazejowski et al. 1 2000), Arbeiter et al.| ( |2002) , 



|Sikora et al] 42009) and |Tavecchio et al. 1(2011) 

In Figure|6]we show the results of our model fitted by eye to the 
simultaneous multi- wavelength spectrum of PKS0227. The model 
fits well to the spectrum across all wavelengths with a transition 
region located at a distance 1.5pc from the central black hole with 
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Figure 4. This figure shows the energy density of the different radiation sources in the rest frame of the jet at different distances along the jet for Fit 1 
to PKS0227. We see that the acceleration of the jet results in a steady increase in the Doppler-boosting of the BLR, dusty torus, NLR, CMB and starlight 
up to the transition region. At different distances along the jet different external photon sources become dominant as we expect. At the transition region 
{x = 1.02 X 10 18 m) the NLR is the dominant energy density resulting in the highest energy gamma-ray emission, however, the CMB becomes the dominant 
energy density at the largest distances along the jet where the majority of inverse-Compton emission is produced. 



Parameter 


Fit 1 


Fit 2 


Fit 3 




1.2 x 10 39 W 


1.5 X 10 39 W 


1.2 X 10 39 W 




8.3 x lO-^W 


8.3 x 10 M W 


8.3 x lO^W 




3.4 x 1O 9 M 


1.5 x 1O 8 M 


2.0 x 1O 9 M 


L 


2 x 10 21 m 


10 21 m 


2 x 10 21 m 


#obs 


0.8° 


0.8° 


1.5° 


^opening 


3° 


3° 


3° 


^-cqui { x — 0) 


0.01 


0.01 


0.01 


M 


3.4 x 1O 9 M 


1.5 x 1O 8 M 


4.5 x 1O 7 M 


a 


1.9 


1.7 


1.5 


-^min 


10.2 McV 


15.3 MeV 


102 MeV 


E 


7.22 GcV 


643 McV 


435 MeV 


70 


4 


4 


4 


7max 


45 


48 


15 


Tmin 


35 


38 


10 


^outcr 


Ikpc 


Ikpc 


Ikpc 



Table 1. This table shows the values of the physical parameters used in the model fits to the multi-wavelength observations of PKS0227 as shown in Figuresp] 
and|6]and[8] 



a jet radius of 9.0 x 10 m. A transition region at this distance 
corresponds to a black hole mass of 1.5 x 1O 8 M if the jet has 
the same geometry as M87 scaled linearly with black hole mass. 
In Table [T] we show the parameters used in the fit. Most of the jet 
parameters are very close to the values found for the fit at large 
distance in the previous section. Again the fit requires a high power 
(Wj = 1.5 x 10 39 W), high bulk Lorentz factor jet (7 max = 48) 



observed close to the jet axis ^obs = 0.8°. The maximum bulk 
Lorentz factor is 48, which decelerates to 38 by x = 10 21 m. 

In Figure [7^ we show the different components of the inverse- 
Compton emission for the fit. We see that the dusty torus photons 
produce the highest energy gamma-rays with the CMB responsi- 
ble for the low energy x-rays as in fit 1. Due to the relatively high 
magnetic field strength in the transition region, 4.81 x 10 J T, the 
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Figure 5. This figure shows the evolution of the electron population contained in a slab of width lm in the plasma rest frame as it travels along the jet in 
Fit 1. At the base of the jet, at x = 0, the electron population given by Equation |25| the electrons suffer large radiative and adiabatic losses in the compact 
accelerating region due to the high magnetic field strength and so the population quickly cools with the highest energy electrons cooling fastest. The population 
cools until it reaches the shock at the base of the conical section at 10 5 iJ s = 1.02 X 10 18 m where additional electrons are accelerated and the plasma is in 
equipartition. The population then proceeds to cool radiatively as it moves away from the shock but with a substantially greater radiative lifetime than at the 
base of the jet. This behaviour is repeated as the electrons are accelerated in the decelerating part of the jet. 



jet produces SSC emission which contributes to the lowest energy 
x-ray observations. In Figure[7]5 we show the emission from differ- 
ent distances along the jet. The synchrotron emission is dominated 
by the transition region as we expect from the discussion in section 
8, with the radio emission coming from larger distances along the 
jet. The high energy inverse-Compton emission is dominated by the 
transition region within the dusty torus, whilst the x-ray emission 
originates from scattering CMB photons at large distances. 

10.3 Fit 3 - A transition region within the BLR 

Finally we investigate modelling the spectrum using a transition 
region within the BLR. An emission region within the BLR has 
been investigated by man y authors including |Sikora et alT] l |1994), 
Bottcher & Bloom|p000t , |Ghisellini & Tavecchio| l |2008) and |Der-| 
mer et al. (2009), and is the most popular location for the gamma- 
ray emission in blazars due to the ability to reproduce short flaring 
timescales. 

In Figure[8]we show a fit to the spectrum of PKS0227 using a 
transition region within the BLR and dusty torus. We show the pa- 
rameters fitted by eye to the observations in Table [T] The transition 
region is located at a distance 0.45pc from the black hole with a 
jet radius of 2.7 x 10 14 m. This corresponds to a black hole mass 
of 3.0 x 10 7 Mq if the jet has the same geometry of M87 scaled 
linearly with black hole mass. We see that whilst the model is able 
to reproduce the high energy gamma-rays and optical synchrotron 
it is unable to fit the radio or x-ray observations. We find that due 
to the short radiative lifetime of electrons in the BLR most of the 



jet power is quickly radiated away leaving insufficient energy in the 
jet at larger distances outside the BLR to reproduce the radio emis- 
sion (see Figure[9j)). This is because the energy in particles is being 
depleted quickly and in order to maintain equipartition the remain- 
ing magnetic energy is constantly converted into replenishing the 
electrons, resulting in most of the jet power being radiated away. 

The high energy inverse-Compton emission is produced by 
scattering BLR photons with a low maximum electron Lorentz fac- 
tor of 851 (see Figure [9^). The high magnetic field strength at the 
transition region due to the small jet radius, 4.59 x 10 _4 T, is re- 
quired to reproduce the high energy synchrotron peak for this low 
maximum electron energy. However, this high magnetic field also 
results in a large amount of SSC emission at x-ray energies which 
in turn requires a very high minimum electron Lorentz factor of 200 
to fit to the x-ray data and not overproduce the observed x-ray emis- 
sion. This very narrow range of electron energies is also difficult to 
understand in the context of shock acceleration where electrons are 
accelerated from low to high energies in a power law. 

The short radiative lifetime within the BLR due to the high 
magnetic field strength and large external photon energy density 
means that any energy contained in high-energy electrons in the jet 
is quickly radiated away. This means that either the jet is still mag- 
netically dominated at these distances and that only occasional tran- 
sient bursts of intense particle acceleration occur within the BLR 
(possibly causing the observed short-timescale gamma-ray flares) 
so that the majority of the jet's energy is carried to larger distances 
as required to reproduce radio observations and by the energy re- 
quirements of radio lobes. 
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Figure 6. This figure shows the results of fitting the model to the SED of PKS0227 using a transition region within the dusty torus. The model fits the 
observations well across all wavelengths including radio observations. We find the model fit requires a high jet power, high Lorentz factor jet observed at a 
small angle to the jet axis. We find that the inverse-Compton peak is due to scattering of Doppler-boosted dusty torus photons at high energies and CMB 
photons at low energies. 
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(a) The fit to PKS0227 showing the different contributions to the inverse- 
Compton emission. 
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(b) The fit to PKS0227 showing the emission from different distances 
along the jet 



Figure 7. This figure shows the different inverse-Compton components in Figure a and the jet emission from different distances in Figure b, using a transition 
region within the dusty torus (Fit 2). Figure a shows that the inverse-Compton emission is produced by scattering dusty torus photons at high energies whilst 
the x-ray emission is due to scattering CMB photons. Figure b shows that the optically thin synchrotron emission is dominated by the transition region, with 
the radio emission coming from larger distances as we expect. The high energy inverse-Compton emission originates close to the transition region within the 
torus whilst the x-ray emission due to scattering CMB photons comes from further along the jet. 



We find that the inferred black hole mass of a transition region 
within the BLR is relatively low and the high peak frequency of the 
corresponding accretion disc radiation makes fitting the spectrum 
difficult. In this fit we use a separate black hole mass for the ac- 
cretion disc Mace and jet geometry M. For the value of M acc we 



use 2 x 10 M© found by |Ghisellini et al.| l |2009| > from fitting an 
accretion disc to the spectrum. 
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Figure 8. This figure shows the results of fitting the model to the SED of PKS0227 using a transition region within the BLR. The model fits the high energy 
synchrotron and gamma-ray observations reasonably well but is unable to reproduce the radio observations due to the short radiative lifetime of electrons within 
the BLR. The fit requires a low maximum electron Lorentz factor to match the gamma-ray peak frequency and a high minimum electron Lorentz factor in order 
to not overproduce the observed x-ray emission with SSC. We find that the inverse-Compton peak is due to scattering of Doppler-boosted BLR photons at high 
energies and SSC photons at low energies. 




(a) The fit to PKS0227 showing the different contributions to the inverse- (b) The fit to PKS0227 showing the emission from different distances 
Compton emission. along the jet 

Figure 9. This figure shows the different inverse-Compton components in Figure a and the jet emission from different distances in Figure b, using a transition 
region within the BLR (Fit 3). Figure a shows that the inverse-Compton emission is produced by scattering BLR photons at high energies whilst the x-ray 
emission is due to SSC. Figure b shows that the optically thin synchrotron emission is dominated by the transition region, with an absence of radio emission 
from larger distances. Both the high and low energy inverse-Compton emission originate from close to the transition region within the BLR due to the high 
magnetic field strength. 



10.4 Constraining the location of the emission region 

We find that both a transition region at large distances and a transi- 
tion region within the dusty torus are able to reproduce the spectrum 
of PKS0227 well. However, we find that a transition region (where 
the jet comes into equipartition in our model and starts emitting) 



within the BLR is unable to reproduce the observations due to the 
short radiative lifetimes of electrons within the BLR. We find that 
the jet must still be magnetically dominated within the BLR since 
the plasma is unable to maintain a state close to equipartition due to 
high radiative energy losses which quickly drain the jet of energy. 
This means it is unable to produce sufficient radio emission at large 
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distances and is also incompatible with the proportion of jet energy 
deposited in radio lobes. It is possible that occasional bursts of par- 
ticle acceleration occur within the BLR (which could be responsible 
for short timescale gamma-ray flares). 

We find that both fits to the spectrum of PKS0227 using a tran- 
sition region at large distances and within the dusty torus require 
large jet power, high bulk Lorentz factor jets observed close to the 
line of sight. This is consistent with Compton-dominant sources be- 
ing FRII type jets. For both these fits we find that the x-ray emission 
is due to inverse-Compton scattering of the CMB at large distances. 
This is the first time that the CMB has been investigated in a model 
for blazar emission so this is an interesting result. 

The black hole mass corresponding to a transition region lo- 
cated at 10 5 i? s as in M87 is 3.4 x 10 9 M o and 1.5 x 10 s M for 
fits 1 a nd 2 respectively. A previous investigation by Ghisellini et al. 
1 2009 ) found a black hole mass of 2 x 10 9 M© by modelling an ac- 
cretion disc spectrum to PKS0227. This black hole mass is close to 
the value of fit 1 which uses a transition region at large distances 
but is substantially larger than fit 2 with a transition region inside 
the dusty torus. 

Whilst our mass estimate of the black hole inferred from fit 1 
is close to that found by Ghisell mi et al.| ([2009 1 the values of our 
jet parameters differ considerably from theirs. This is due to funda- 
mentally different approaches to modelling the jet. In their model 
they consider a single compact blob which produces the observed 
inverse-Compton emission from scattering BLR photons whilst the 
optical observations are matched using thermal accretion disc, torus 
and corona components. Since they only consider a single emitting 
region they do not attempt to reproduce the radio data. In this model 
they find a bulk Lorentz factor of 14 substantially smaller than in 
our jet, a viewing angle is 3° slightly larger than in our fit and a total 
jet power is higher by over an order of mag nitude 2.38 x 10 40 W, 
however, their model includes one proton per electron. These very 
different parameters are due to fundamentally different model as- 
sumptions. It is worth pointing out that no attempt is made to fit the 
model to the observed radio data and in addition the synchrotron 
radiation also does not fit the optical data since this is fitted by the 
thermal components, so the non-thermal synchrotron emission is 
not well constrained. 

The geometry of M87 suggests that the transition region where 
the jet plasma approaches equipartition is at relatively large dis- 
tances from the black hole, > lOpc for black holes of mass > 
10 9 Mq. This puts the transition region outside the BLR and out- 
side or near to the outer edge of the dusty torus so we do not expect 
photons from these regions to dominate the external radiation field 
at these distances (see Figure[4|. We find that this is consistent with 
our model of the observed spectrum of PKS0227 which is well fit- 
ted by a transition region at a distance 34pc from the black hole 
where the inverse-Compton emission is due primarily to scattering 
of CMB photons. The transition region may not occur at 10° R s for 
every jet as it could depend on the jet power, mass loading etc. If 
this were the case we might expect a jet with a high bulk Lorentz 
factor to have a transition region at a larger distance than the com- 
paratively low power jet in M87, in order to accelerate to a higher 
bulk Lorentz factor. This is not the case for fit 2 where the transition 
region would have to be located at 7500i? s to be compatible with 
the black hole mass estimate of Ghisellini et al. ( 2009 1. For these 
reasons we favour a scenario in which the transition region occurs 
at large distances outside the BLR and dusty torus. 

From the Figures[2]and[6]we see that one of the distinguishing 
features between a transition region at large distances or within the 
dusty torus is the location of the optically thick to thin synchrotron 



break. Since the optically thin synchrotron emission is dominated 
by the transition region (see Figures [5J) and[7J)) the location of the 
transition region largely determines the frequency at which the opti- 
cally thick to thin break occurs. We see that the frequency at which 
the flat self-absorbed radio spectrum becomes optically thin and 
changes gradient increases for a transition region at a smaller dis- 
tance along the jet as we expect due to the increased magnetic field 
and smaller jet radius. This feature could potentially be used to dis- 
tinguish between a transition region at large distances or within the 
dusty torus. We wish to investigate whether this can be used to con- 
strain the emission region in the near future. 

10.5 An adiabatic or ballistic conical jet? 

In this paper we have assumed that the parabolic accelerating re- 
gion of the jet expands adiabatically whilst the conical region is 
ballistic. In Paper I we proposed that a continuous conical jet with 
a constant opening angle and bulk Lorentz factor would not suffer 
adiabatic losses because it does no work expanding its environment 
(since the volume occupied by the jet is constant with time in the 
lab frame). In an idealised homogeneous jet the velocity profile in 
the direction perpendicular to the jet axis scales linearly with expan- 
sion and so small sections of the jet maintain a constant perpendic- 
ular velocity and are effectively ballistic. Whilst these arguments 
are based on an idealised situation, a conical jet with an opening 
angle Opening 7buik > 1 is not causally connected laterally and an- 
alytic work by Zakamska et al. (2008 1 found that far away from the 
black hole a conical jet will always become causally disconnected 
and this seems to be the case for jet parameters inferred from ob- 
servations of GRBs [Kohler & Begelman (2012). Work reproducing 
the large scale observed emission and polarisation of FRI jets using 
MHD models by |Laing & Bridle| f2004[ > found that jets were bet- 
ter described by non-adiabatic jets as opposed to adiabatic jets. It 
is therefore worth testing the effect of including adiabatic losses in 
the conical section of our model. 

In Figure [To] we show the results of including adiabatic losses 
on fits 1 and 2. We find that in fit 1 adiabatic losses cool the high 
energy electrons more quickly than scattering the CMB and so the 
inverse-Compton emission is vastly reduced. In the case of adia- 
batic losses we find it difficult to fit the spectrum of PKS0227 with 
a transition region at large distances without using very large bulk 
Lorentz factors. In the case of fit 2 we also find that adiabatic losses 
substantially decrease the inverse-Compton emission from the jet 
especially the x-ray emission which was due to scattering CMB 
photons. We also find that adiabatic losses result in a decrease in the 
radio emission at large distances due to the adiabatic energy losses 
along the jet. We find it difficult to fit the spectrum of PKS0227 at 
radio and x-ray energies using a transition region within the dusty 
torus including adiabatic losses. 

For these reasons we favour a ballistic conical section of the 
jet. This could be a result of the assumptions we have made when 
fitting the spectrum with our jet model. We have assumed that the 
jet comes into equipartition at the transition from parabolic to con- 
ical and the jet maintains equipartition outwards from this region. 
This is motivated by simulations which indicate the jet transitions 
to conical when the plasma approaches equipartition and from ob- 
servations which seem to indicate that the jet plasma is observed 
to be close to equipartition with fairly continuous optical emission 
requiring continuous particle acceleration along the jet. In this pa- 
per we have assumed a simple scenario in which the jet suddenly 
transitions to equipartition, however, it is likely that some particle 
acceleration and emission comes from the magnetically dominated 
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(a) Fit 1 to PKS0227 with and without adiabatic losses in the conical part (b) Fit 2 to PKS0227 with and without adiabatic losses in the conical part 
of the jet. of the jet. 

Figure 10. This figure shows fits 1 and 2 to the spectrum of PKS0227 with and without adiabatic losses in the conical section of the jet. The Figure shows that 
the result of adiabatic losses on the electron population is to severely reduce the amount of inverse-Compton radiation, especially from scattering CMB photons 
at large distances. The radio emission is reduced as the jet loses additional energy to adiabatic losses at large distances where the radio emission originates. We 
find it difficult to reproduce the spectrum using a transition region at large distances with adiabatic losses without using a large bulk Lorentz factor. We also 
find it difficult to reproduce the spectrum using a transition region within the dusty torus due to the reduction in x-ray and radio emission. 



region of the jet and this could affect our conclusions. We intend to 
investigate different mechanisms for particle acceleration along the 
jet in the near future. 

We have found that the x-ray emission is well fitted by inverse- 
Compton scattering of CMB photons at large distances along the jet. 
This x-ray emission is severely reduced if the conical section is adi- 
abatic and becomes difficult to reproduce by non-thermal emission 
from the jet. The x-ray emission may also be due to the accretion 
disc corona so an inability to reproduce the x-ray emission with 
the jet does not necessarily exclude an adiabatic conical jet. We 
are keen to investigate whether this result holds for other Compton- 
dominant jets in the next paper in the series. 

10.6 Discussion 

In this work we have attempted to model the quiescent state of the 
blazar PKS0227 since this constrains the average jet parameters. 
Short term variability is observed on time-scales as short as hours 
or minutes for the most highly variable blazars (Aharonian et al. 
( |2007| >). The quiescent emission from our jet model comes from 
large distances with long electron cooling times ( ~years for fit 
1 and ~days for fit 2), so is not easily compatible with this short 
timescale variability. In this series of papers we wish to initially 
investigate the average quiescent emission from different elements 
of the blazar population in order to check that our model is able to 
reproduce their spectra and to find any trends in the quiescent jet pa- 
rameters. Once we are confident that our model is able to reproduce 
the quiescent spectra of both BL Lacs and FSRQs with reasonable 
jet parameters, we will investigate whether we can reproduce flar- 
ing spectra by allowing flares to occur within the parabolic region of 
the jet. The parabolic region of the jet is more compact with shorter 
cooling timescales due to the higher magnetic field strength and en- 
ergy density of external photons, so this region is more likely to be 
responsible for short timescale flaring. 

It is important to emphasise that whilst our model is more com- 
plicated than a single zone model it does not possess many more 
free parameters. We have attempted to minimise the number of free 
parameters in our model by making reasonable, simplifying phys- 



ical assumptions and using the results of previous work observing 
and modelling the external radiation fields. The CMB is well con- 
strained by observations and possesses no free parameters, how- 
ever, the emission and distribution of the NLR is less well known. 
We have fixed most of the free parameters using values from previ- 
ous observations and modelling of the BLR, dusty torus and NLR 
and we have left only the outer radius a; ou tor of the NLR as a free 
parameter since this seems to vary the most significantly. In these 
fits we have not injected an arbitrary electron population into ar- 
bitrary sections along the jet, we have injected a small amount of 
the initial jet power (1%) into electrons at the very base of the jet 
(to illustrate the resulting spectrum from electrons at the base of 
the jet whilst not affecting the fit to observed data, see Figure [3]}). 
The acceleration of electrons in the decelerating section of the jet 
is determined entirely by Equation [T7] in terms the other jet param- 
eters in Table [T] and assuming equipartition in the conical section. 
Along with the fixed geometry this leaves our model with twelve 
free parameters which we vary when fitting the spectrum (Wj, L, 

Lacc> ^obst ^opening* -£/rnin; -E^max; 'Y max ' ^fmin and Xouter)- 

This number of free parameters is comparable to single-zone mod- 
els with broken power-law electron distributions. 

We have tested the effect of changing the functional form of 
the deceleration of the bulk Lorentz factor along the jet (Equa- 
tion |5j. In our model the bulk kinetic energy of the jet is spread 
evenly on a logscale along the jet from deceleration. We find that 
strongly decelerating the jet over a particular distance scale effec- 
tively boosts the relative amount of emission from that section of 
the jet. This means that if the jet decelerates strongly over a scale 
~ 10 5 i? s then the amount of optical synchrotron and SSC emis- 
sion is boosted relative to the radio emission and inverse-Compton 
scattering of CMB photons. Similarly if the jet is chosen to decel- 
erate strongly at larger scales the amount of low frequency radio 
and inverse-Compton emission increases relative to the optical syn- 
chrotron and SSC emission. We find that strong deceleration on an 
intermediate length scale leads to a bump in the characteristic radio 
frequencies of that section. This deviation from a flat radio spec- 
trum is not something seen in radio observations of blazars and so 
we disfavour strong deceleration of the jet on an intermediate length 
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scale 10 17 — 10 19 m. We find that deceleration which decreases as 
the logarithm of the jet distance as we have assumed (Equation [5} 
distributes power more evenly along the jet resulting in a radio spec- 
trum which is smooth and close to flat consistent with observations. 

We find that electrons injected into the parabolic section of the 
jet quickly radiate all their energy away due to the large magnetic 
fields and high energy density in the seed photon field (see Fig- 
ure^). We also find that the emission of these electrons tends to 
be dominated by the emission from the conical section of the jet, 
since the inner regions of the accelerating parabolic section have 
a lower bulk Lorentz factor and so are not so strongly Doppler- 
boosted. Since all the energy injected into electrons in the inner 
parts of the parabolic region is quickly radiated away and does not 
contribute as significantly to the observed spectrum it is difficult to 
constrain the amount of energy injected into accelerating electrons 
in this region. Since we expect the jet to be magnetically dominated 
in the parabolic region we have assumed that a small percentage of 
the magnetic energy goes into shock acceleration of the electrons at 
the base of the parabolic section, this means that our estimated jet 
power is a lower bound because the total jet power would increase if 
a larger proportion of magnetic energy were converted into particle 
energy in the inner parabolic region, without a significant increase 
in observed emission. 

We find that the low frequency radio emission of the jet is dom- 
inated by the conical sections with the lowest bulk Lorentz factors 
and largest radii (see Figure[3j)). This is partly because the regions 
with high bulk Lorentz factors are strongly Doppler-boosted so the 
synchrotron opacity in these sections is effectively the synchrotron 
opacity for a much lower energy radio photon than that observed. 
The synchrotron opacity increases steeply at low frequencies so the 
radio emission from regions with high bulk Lorentz factors is sup- 
pressed and the strongly decelerated regions with larger radii dom- 
inate the radio. We see that the model is able to fit to the nearly 
flat radio spectrum of PKS0227 very well. Including the effects of 
deceleration and the conversion of bulk kinetic energy into particle 
acceleration means that the electron energy losses are replaced as 
they travel along the decelerating conical part of the jet as expected 
from observations of optical synchrotron emission far from the base 
of a jet ( |Meisenheimer et al.|{1997| > and |Jester et al.H200l) ). This 
mechanism to rejuvenate the high energy electron population has 
been investigated a number of times previously, however, this is the 
first investigation that has included the effects of deceleration in a 
continuous conical jet which conserves energy-momentum and par- 
ticle number. 

Most previous investigations have neglected radio observa- 
tions and concentrated on one or two-zone models which are dis- 
connected from the large scale radio emission of the jet. We have 
demonstrated that our rigorous, physically motivated, inhomoge- 
neous jet model fits well to the simultaneous multiwavelength spec- 
trum of PKS0227 across all wavelengths including radio observa- 
tions, which are usually ignored. This allows us to better constrain 
the regions where emission may occur in the jet by reproducing the 
radio and large scale emission of the jet whilst conserving energy 
along the jet. We find the results of this initial application of our 
model to a blazar very promising and we intend to use this model 
to conduct a more detailed study of a sample of Compton-dominant 
blazars in the near future to see if our results hold for other blazars 
and to try to distinguish between a transition region at large dis- 
tances or within the dusty torus. 



11 CONCLUSION 

In this paper we have developed a realistic jet model which in- 
cludes an accelerating, magnetically dominated, parabolic base 
which transitions to a slowly decelerating conical jet with a jet ge- 
ometry set by recent radio observations of M87 and consistent with 
simulations and theory. This is the first time that an accelerating 
parabolic region transitioning to a decelerating conical region has 
been treated in a jet model for the spectra of a blazar. We calculate 
the observed line of sight synchrotron opacity to each section along 
the jet by Lorentz transforming into the rest frame of the plasma. We 
calculate the inverse-Compton emission from accretion disc pho- 
tons directly from the disc, BLR, dusty torus, NLR, starlight and 
CMB photons by Lorentz transforming the radiation into the plasma 
rest frame. Our model conserves relativistic energy-momentum and 
particle number along the jet and we take into account radiative and 
adiabatic losses on the electron population. 

We assume that the jet is initially magnetically dominated until 
it transitions from parabolic to conical at 10 5 7? s , where the jet stops 
accelerating and the plasma approaches equipartition. Electrons are 
accelerated at this transition region and electron energy losses are 
replenished along the jet by the conversion of bulk kinetic energy 
into particle energy as the jet decelerates so that the plasma is kept 
in equipartition throughout the conical section of the jet. We use the 
model with a geometry and acceleration set by VLBI observations 
of M87 and simulations to fit to the simultaneous multi-wavelength 
spectrum of the Compton-dominant FSRQ PKS0227. 

We investigate a transition region located outside the dusty 
torus, within the dusty torus and within the BLR. We find that our 
model fits very well to the spectrum across all wavelengths includ- 
ing radio observations with a transition region at either large dis- 
tances from the central black hole, or within the dusty torus. We 
find that due to the short radiative lifetime we are unable to re- 
produce the radio observations using a transition region within the 
BLR. We find that within the BLR the jet must still be magnetically 
dominated in order for the jet to transport sufficient energy to larger 
distances, however, this does not exclude occasional transient bursts 
of intense particle acceleration. 

We find that our fits using transition regions outside the BLR 
both require high power (~ 10 39 W), high bulk Lorentz factor jets 
(> 40) viewed at a small angle to the jet (~ 1°). We find an in- 
ferred black hole mass of 3.4 x 1O 9 M0 for the transition region 
outside the dusty torus and 1.5 x 10 s M© for a transition region 
within the dusty torus using observations of the geometry of M87 
scaled linearly with black hole mass. The inferred black hole mass 
of the transition region at large distances, which completely deter- 
mines the i nner geometry of the j et, is close to the black hole mass 
found by |Ghisellini et al.| ([2009) (2 x 10 9 A/ Q ). We find that the 
optically-thick to thin synchrotron break depends on the location 
of the transition region and can potentially be used to distinguish 
between the two emission locations. 

Interestingly, we find that including the redshift dependence of 
the CMB and by boosting into the plasma rest frame the Compton- 
dominance of the transition region at large distances 34pc is well 
fitted by inverse-Compton scattering of CMB photons. This could 
explain why Compton-dominant blazars are powerful, high bulk 
Lorentz factor blazars which are preferentially at high redshift due 
to the temperature dependence of the CMB oc (1 + z). We find that 
the x-ray emission for both transition region fits outside the BLR 
can be fitted well by scattering of CMB photons on large scales if 
the conical section of the jet is ballistic. 

We find these initial results very promising and we intend to 
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use this model as the basis for a more detailed investigation into the 
physics behind different elements of the blazar population in the 
near future. 
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